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Abstract 



We study the analytic structure of the S-matrix which is obtained from 
the reduced Wheeler-DeWitt wave function describing spherically symmetric 
gravitational collapse of massless scalar fields. The complex simple poles in 
the S-matrix lead to the wave functions that satisfy the same boundary con- 
dition as quasi-normal modes of a black hole, and correspond to the bounded 
states of the Euclidean Wheeler-DeWitt equation. These wave function are 
interpreted as quantum instantons. 
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In the previous work [[!]] we studied quantum mechanically the self-similar black hole 
formation by collapsing scalar fields and found the wave functions that give the correct 
semi-classical limit. The reduced Wheeler-DeWitt equation for gravity belongs to an exactly 
solvable Calogero type system with an inverted potential whose attractive inverse square and 
repulsive square potential terms give rise to a potential barrier. The boundary condition for 
black hole formation was that the wave function has both the incoming and the outgoing 
flux at spatial infinity and only the incoming flux toward the black hole singularity. 

Of particular interest is the subcritical case, in which a black hole can be formed through 
quantum tunneling. Due to the time reversal symmetry, however, the subcritical wave 
function may be given an interpretation of the reversal process of black hole formation, that 
is, the decay of the black hole |§. Then the wave function for black hole decay should have 
a purely outgoing flux. This wave function is somehow reminiscent of the gravitational wave 
from a perturbation of a black hole || . Moreover, for a certain discrete spectrum of complex 
frequencies there occur quasi-normal modes that have both purely outgoing modes at spatial 
infinity and purely incoming ones at the horizon of the black hole |4]]. 

In this paper we study the pole structure of the S-matrix which is obtained from the 
wave function for black hole formation. The boundary condition that the wave function 
should have a purely outgoing flux at spatial infinity and a purely incoming one at the 
classical apparent horizon leads to a discrete spectrum of complex parameters (cq). It is 
further shown that this wave function can be obtained through the analytical continuation 
of a bounded state of the corresponding Euclidean Wheeler-DeWitt equation. Just as quasi- 
normal modes of perturbations of a black hole can be interpreted as instantons ||, these 
exact wave functions of the quantum theory for black hole decay may be interpreted as 
quantum instantons 

The spherically symmetric geometry minimally coupled to a massless scalar field is de- 
scribed by the reduced action in (1 + l)-dimensional spacetime of which the Hilbert-Einstein 
action is 
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The reduced action is 

Ssph ~ 

where jab is the (1 + l)-dimensional metric. The spherical spacetime metric is 

ds 2 = —2du dv + r 2 dQl, 
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where dVt\ is the usual spherical part of the metric, and u and v are null coordinates. The 
self-similarity condition is imposed such that 



r = y/-uv y(z), <f> = <f>(z), 
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where z = +v/(—u) = e 2r , y and <fi depend only on z. We introduce another coordinates 



u = —toe , v — toe 
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to rewrite the metric as 



ds 2 
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where N(t) is the lapse function of the ADM formulation. 

The classical solutions of the field equations were obtained by Roberts 0, and studied 
in connection with gravitational collapse by others 0. Classically black hole formation is 
only allowed in the supercritical case (c > 1), but even in the subcritical situation there 
are quantum mechanical tunneling processes to form a black hole of which the probability 
is semiclassically calculated 

In our previous work |jj we quantized the system canonically with the ADM formulation 
to obtain the Wheeler-DeWitt equation for the quantum black hole formation 
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where K/% = (m 2 /h 2 )(u 2 /2) plays the role of a cut-off parameter of the model, and we use 
a unit system c = 1. The wave function can be factorized to the scalar and gravitational 
parts, 
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Here the scalar field part is chosen to yield the classical momentum = hKy 2 <p/m 2 pN = 
±Kcq, where cq is the dimensionless parameter determining the supercritical (cq > 1), the 
critical (co = 1), and the subcritical (1 > cq > 0) collapse. 

Now the gravitational field equation of the Wheeler-DeWitt equation takes the form of 
a Schrodinger equation 



-h 2 d 2 K ( 

~2Kdf + Y{ 2 



y I 



1>(y) = o. 



The solution describing black hole formation was obtained in Ref. [Q: 
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where M is the confluent hypergeometric function and 

a=--^-(Q + K), b=l 
with 
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Using the asymptotic form M at spatial infinity 
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we obtain the S-matrix component describing the reflection rate 

s _ T(b - a) {iK/h) 2a - b 
~ T(a) e ina 

From the S-matrix follows the transmission rate for black hole formation 

Jtrans 
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is used. Equation (\T^) gives the probability of black hole 
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formation for the supercritical, critical, and subcritical Co-values. 

We now consider the analytic structure of the S-matrix: it is an analytic function of Q 
and K with simple poles which can be explicitly shown as 
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The poles reside in the unphysical region of the parameter space of Q and K: 

Q-K = -ih(2N + l), (N = 0,1,2, -••). 
It should be remarked that these poles make the first term of Eq. ( |i~3| ) vanish since 

b-a = ~(Q-K) + ± = -N, (N = 0, 1, 2, ■ ■ ■). (18) 
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The second term of Eq. ( p~3| ) has a purely outgoing flux at spatial infinity. The wave function 
near the apparent horizon, which can be obtained by the steepest descent method in the 
Appendix of Ref. M and by taking the large (iT/7i)-limit, leads to the flux 
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where A(y) denotes an amplitude, a real function, and Cq 



i(^/i^)(2A^ + 1) from Eq. 



(|l7l) and c = (eg — Ti 2 / AK 2 ) 1 ! 2 — i{%/K). At the apparent horizon y^u = Co/a/2, the wave 
function has an incoming flux. Therefore, the poles are the outcome of the same boundary 
condition used to find quasi-normal modes of a black hole [|J. Note that the wave function 
(|j~0|) has also the purely incoming flux toward the black hole singularity at y — 0. 

A few comments are in order. First, for physical processes of gravitational collapse there 
can not be poles because K and Co are real-valued. In ordinary quantum mechanics, the poles 
of S-matrix occur at the bound states [ 1Q| , and in relativistic scattering at the resonances or 
the Regge poles [|ll|]. Our case is analogous to a meta-stable quantum mechanical system of 
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which poles are identified with quasi-stationary states that describe the decay of a particle 
through a potential barrier. Second, we calculated quantum decay rate of a black hole 
as a reversed process of gravitational collapse through a barrier by quantum tunneling. 
This quantum decay process, first studied by Tomimatsu @, is a distinctively different 
decay channel from the Hawking radiation process. It will be interesting to investigate both 
processes present in one model. Third, it should be pointed out that our discussion based 
upon the similarity of the boundary condition on the wave function with quasi-normal modes 
seems to have no deeper physical connection more than analogy because our model works 
only for a dynamical stage of gravitational collapse and its reversed process, rather than the 
quasi-stationary stage at late times. 

Recalling that the poles in Eqs. ( |16|) or (|i~8l) result from the potential barrier and that the 
exponential behavior of the Wheeler-DeWitt equation under a potential barrier describes 
a Euclidean geometry, we turn to the Euclidean theory of gravitational collapse. In the 
Euclidean theory the Wheeler-DeWitt equation has oscillatory wave functions and a well- 
defined semiclassical limit even under the potential barrier of the Lorentzian theory [fl2 
The Euclidean geometry with the metric 
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leads to the Wheeler-DeWitt equation 
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the Lorentzian and Euclidean geometries [12], the wave function has the form 
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The Wheeler-DeWitt equation reduces to the gravitational field equation 

-h 2 d 2 K( 2 c|_ N 
2K dy 2 + 2 [ V + y 2 , 
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Notice that this is a variant of Calogero models with the Calogero-Moser Hamiltonian |L3 



but the energy eigenvalue is fixed, and only a quantized c is allowed. Since Eq. ( |23"D can 
also be obtained from the Lorentzian equation (H) by letting 

K = iK E} (24) 

one may obtain, through the analytical continuation of Eq. (0), the solution to Eq. ( p3|) 
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where 



a E = 2 + 2ft{QE + K E ) 
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with 

Qe= Klcl + ^) . (27) 



The asymptotic form of Eq. flzq) leads to the bounded states only when 

b E -a E = -N, (N = 0,1,2, ■■■), (28) 

that is, the condition is satisfied 

Q E -K E = -h{2N +1). (29) 

The condition (|29|) is identical to the pole position of the S-matrix with K = iK E given in 
Eqs. (|T7p and flE|). 

A few remarks are in order. First, the quantum solution (|25|) is analogous to an instanton 
in the sense it is a solution in the Euclidean sector, but is not in the strict sense because 
the Wheeler-DeWitt equation is already a quantum equation, not a classical one. The 
semiclassical result from the Bohr-Sommerfeld quantization rule 

Ix (1_C ° )=7r ( iV + ^' ( iV = ' 1 ' 2 '---), (30) 



is the large (K E /h)-\imit of the exact result (p29|) . The instantons, the left hand side of 
Eq. (|3Tj|), provides semiclassically the probability of tunneling process ||. Second, the 
correspondence between the poles and the Euclidean polynomial solutions breaks down for 
large N. While the poles contribute for all N without limit, the normalizable Euclidean 
solutions exist only for N < K E /2h. The polynomial solutions for large N are well defined, 
but are not normalizable. We have not yet understood these nonnormalizable solutions. 
Finally, we consider the classical field equations corresponding to the poles of the S-matrix. 
In the Lorentzian geometry the relevant equations are 

dr v 



K 2 [ —2 + y 2 + -^r) , (32) 





where Co — 1 — ifi(2N + 1)/K, for large K. The complex cq implies complex ^ and 
which may be imagined as a bound state like complex momentum in quantum mechanics 
and requires a complex spacetime metric. In the Euclidean geometry (K = iK E ) these 
classical equations are the same as those equations with quantized c in the tunneling region 
in Ref. HI. 
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